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Abstract 

This  study  investigated  the  atmospherie  reentry  of  the  Spaee  Maneuver  Vehiele  from 
low-earth  orbit  using  an  entry  guidanee  eoneept  similar  to  the  Spaee  Shuttle.  The  Spaee 
Maneuver  Vehiele  was  modeled  as  a  point  mass  with  aerodynamie  properties  as 
determined  using  Newtonian  impaet  theory.  For  the  rarefied-flow  transition  regime 
bridging  formulae  are  used  to  eapture  the  effeets  of  both  hypersonie  eontinuum  flow  and 
free  moleeular  flow.  Constraints  to  the  reentry  are  developed  and  a  reentry  eorridor  is 
defined  in  the  drag-veloeity  plane.  Bank  angle  modulation  is  established  as  the  primary 
means  for  eontrolling  drag  and  range.  The  guidanee  eoneept  is  applied  to  both  a  high 
inelination  orbit  and  a  low  inclination  orbit  with  crossrange  requirement.  Monte  Carlo 
error  analysis  validates  the  ability  of  the  control  algorithm  to  guide  reentry  in  the 
presence  of  initial  state  errors,  as  well  as  atmospheric  variations. 
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INVESTIGATION  OF  ATMOSPHERIC  REENTRY 
FOR  THE  SPACE  MANEUVER  VEHICLE 


I.  Introduction 


Background 

Over  the  last  decade,  the  importance  of  space  assets  to  the  warfighter  has  been 
demonstrated  through  various  regional  conflicts.  The  ability  to  communicate  with  troops 
in  the  field,  predict  weather  conditions,  and  monitor  the  enemy’s  movements  from  space 
all  combine  to  give  theater  commanders  an  unprecedented  battlefield  awareness.  Over 
the  next  few  decades  control  of  the  high  ground  of  space  will  play  an  even  more  crucial 
role  in  regional,  as  well  as  global,  conflicts. 

At  the  current  time,  launching  payloads  into  space  is  a  lengthy  process  sometimes 
taking  years  to  complete.  Additionally,  the  combined  cost  of  the  satellite  and  launch 
vehicle  is  typically  in  the  hundreds  of  millions  of  dollars  range.  Not  surprisingly, 
“routine,  reliable,  and  lower-cost  space  operations  are  goals  of  U.S  military  planners” 
(Shaw,  2000:1).  One  concept  from  the  Air  Force  Space  Command  (AFSPC)  Strategic 
Master  Plan  proposed  to  meet  these  goals  is  the  Military  Spaceplane  (MSP)  (Shaw, 
2000:1).  The  MSP  consists  of  four  primary  components  (Shaw,  2000:1): 

•  Space  Operations  Vehicle  (SOV):  a  reusable  launch  element  to  reduce  the 
cost  of,  improve  the  flexibility  of,  and  increase  the  responsiveness  of 
earth-to-orbit  operations; 
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•  Space  Maneuver  Vehicle  (SMV):  an  uncrewed,  reusable  satellite  bus  and 
upper  stage  with  significant  maneuvering  capability; 

•  Modular  Insertion  Stage  (MIS):  a  low  cost  expendable  upper  stage  to 
reduce  the  cost  and  improve  the  operability  of  Space  Support  missions; 

•  Common  Aero  Vehicle  (CAV):  an  aeroshell  designed  to  deliver  material 
(e.g.  munitions,  UAVs,  or  critical  supplies)  through  space  directly  to  a 
theater. 

For  this  study,  the  SMV  component  of  the  MSP  is  of  primary  concern.  The  SMV 
is  envisioned  to  have  between  7000  and  12000  fps  of  A  V— approximately  100  times  the 
maneuvering  capability  of  current  satellites  (Shaw,  2000:6).  The  maneuvering  capability 
combined  with  the  intended  mission  duration  of  less  than  12  months  (on  the  order  of 
weeks  for  tactical  missions)  allows  the  SMV  to  perform  significant  orbital  maneuvers  to 
accomplish  the  intended  mission.  For  reference,  in  a  medium  altitude  earth  orbit  a  AV  of 
10000  fps  will  allow  for  a  20°-plane  change  (Shaw,  2000:7).  In  addition  to  on  orbit 
maneuverability,  the  SMV  is  intended  to  be  responsive  to  the  needs  of  the  warfighter. 

The  goal  is  for  “warfighers  [to]  have  tailored  satellite-based  services  soon  after  they 
requested  them”  (Shaw,  2000:7). 

A  Military  Utility  Analysis  (MU A)  conducted  during  1998  and  1999  showed  that 
the  capabilities  of  the  SMV  would  most  likely  be  used  for  force  enhancement,  space 
control,  and  space  test  in  the  near  future  (Shaw,  2000:9-10).  Comparison  of  the  SMV 
capabilities  against  the  Air  Force  Space  Command  Strategic  Master  Plan  (SMP)  for  2000 
shows  the  SMV  will  have  a  “direct  and  substantial”  effect  on  4  of  the  top  10  priorities  in 
the  near-term,  and  6  of  the  top  10  priorities  in  the  mid-term  (Shaw,  2000:1 1).  According 
to  the  Air  Force  Research  Lab’s  Space  Vehicles  Directorate  (AFRL/SV)  the  SMV  could 
perform  the  following  classes  of  missions  (AFRL  2004): 
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•  Intelligence,  surveillance  and  reconnaissance  of  ground  targets  (with 
either  integrated  or  deployable  ISR  satellite) 

•  Deployment  and  recovery  of  microsatellites  (e.g.  Space  Control  Satellites) 

•  Rapid  constellation  replenishment 

Currently  there  are  two  technology  programs  being  worked  to  demonstrate  the 
feasibility  of  the  SMV  concept.  These  two  programs  are  the  Boeing  and  US  Air  Force  X- 
40A  and  the  NASA  and  Boeing  X-37  (Shaw,  2000;  17).  The  Air  Force  is  also 
contributing  funds  to  the  X-37  program  to  demonstrate  technologies  for  future  reusable 
military  spacecraft  (Shaw,  2000:17). 

The  Air  Force  X-40A  (22  feet  long  and  weighs  about  2600  lb)  is  a  test  vehicle 
designed  to  evaluate  the  low  speed  atmospheric  flight  dynamics  of  the  X-37  (Dry den 
2004).  The  X-40A,  shown  if  Figure  1,  successfully  completed  its  series  of  seven  test 


X-4M  ffM  rigM  M 


Figure  1;  X-40A  Space  Maneuver  Vehicle  (Dryden  2004) 
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flights  at  NASA  Dryden  Research  Center  where  it  autonomously  acquired  the  runway 
and  executed  an  aircraft  landing  after  being  dropped  from  an  Army  Chinook  helicopter  at 
15,000  ft  (X-37  2004).  The  X-40A  is  often  referred  to  as  the  X-40A  Space  Maneuver 
Vehicle,  however,  it  is  a  technology  demonstrator  for  the  Space  Maneuver  Vehicle 
concept.  For  the  purposes  of  this  study,  the  fully  operational  SMV  will  be  considered, 
not  the  X-40A. 

The  X-37  (shown  in  Figure  2)  is  being  designed  to  operate  on  orbit  as  well  as 
during  the  reentry  phase  of  flight.  The  X-37  is  capable  of  fitting  into  the  Space  Shuttle 
payload  bay  for  launch  into  orbit,  or  it  can  be  carried  to  orbit  by  an  expendable  launch 
vehicle  (Marshall  2004).  Technologies  that  will  be  demonstrated  by  the  X-37  include: 
advanced  airframe  design,  avionics  technologies,  and  advanced  thermal  protection 
systems  (Marshall  2004). 


X-37  RsusablB  Spacspldne 

Figure  2:  X-37  Reusable  Spaceplane  (Marshall  2004) 
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The  goal  of  the  current  study  is  to  develop  an  entry  guidance  algorithm  capable  of 
guiding  an  operational  Space  Maneuver  Vehicle  through  atmospheric  reentry  to  a  pre¬ 
determined  termination  point  in  the  vicinity  of  the  runway.  The  currently  operational 
Space  Transportation  Systems  (STS)  Space  Shuttle  Orbiters  are  an  obvious  starting  point 
for  developing  the  guidance  algorithm.  The  Space  Shuttle  has  been  operational  since 
1981  and  has  demonstrated  the  capability  of  its  reentry  guidance  system  on  over  100 
landings.  Therefore,  this  study  will  focus  on  applying  the  concepts  used  in  Space  Shuttle 
entry  guidance  to  the  Space  Maneuver  Vehicle. 

Space  Shuttle  entry  guidance  is  “designed  on  the  principal  of  analytically  defining 
a  desired  drag  acceleration  profile  and  commanding  the  orbiter  to  achieve  the  desired 
profile”  (Harpold  &  Graves,  1978:103).  Thus,  the  guidance  approach  employed  is  to 
establish  the  reentry  corridor  in  the  drag-velocity  (D-V)  plane  by  defining  the  constraints, 
define  a  desired  drag  acceleration  profile  in  the  D-V  plane,  and  then  command  the 
vehicle  to  fly  on  the  defined  profile.  The  primary  control  variable  used  during  the  entry 
phase  is  the  bank  angle,  (j).  The  bank  angle  is  modulated  to  control  both  the  drag 
acceleration  level  and  the  range  to  the  termination  point. 

In  order  to  apply  the  shuttle  entry  guidance  concept  to  the  SMV,  the  reentry 
problem  must  first  be  setup.  This  is  the  focus  of  Chapter  2.  In  Chapter  2,  the  equations 
of  motion  governing  a  point-mass  vehicle  reentering  the  Earth’s  atmosphere  are 
developed.  Next,  since  the  SMV  exists  in  concept  only,  its  aerodynamic  properties  are 
derived. 
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Having  defined  the  problem  and  the  aerodynamie  properties  of  the  vehiele  in 
Chapter  2,  Chapter  3  develops  the  guidanee  algorithm.  First  the  initial  and  final 
eonditions  for  the  reentry  trajectory  are  developed.  Next,  the  entry  corridor  is  defined  by 
deriving  the  constraints  in  the  D-V  plane.  A  reference  profile  in  the  D-V  plane  is  then 
developed.  With  the  reference  profile  within  the  entry  corridor  established,  the  control 
law  and  guidance  algorithm  which  will  guide  the  vehicle  during  its  descent  are 
developed. 

Chapter  4  examines  two  scenarios  of  the  SMV  reentering  the  Earth’s  atmosphere 
with  the  Shuttle  entry  guidance  concept  applied.  Results  are  presented  which  show  the 
entry  profile  within  the  D-V  plane,  as  well  as  the  bank  angle  and  angle  of  attack  profiles. 
Additionally,  the  ground  track  of  the  reentry  path  and  a  plot  of  azimuth  angle  show  the 
effects  of  the  roll-reversal  commands  within  the  guidance  algorithm.  Finally,  a  100-run 
Monte  Carlo  error  analysis  simulation  is  run  to  examine  the  performance  of  the 
developed  guidance  system  in  the  presence  of  errors  in  the  initial  state,  aerodynamics, 
and  atmosphere. 

Chapter  5  then  summarizes  the  results  of  applying  the  concepts  of  Shuttle  entry 
guidance  to  the  Space  Maneuver  Vehicle.  The  results  are  interpreted  to  show  the  success 
of  the  entry  guidance  algorithm.  Areas  of  improvement  within  the  design  are  noted,  and 
recommendations  are  made  for  future  areas  of  study. 

As  a  final  note,  the  author  wishes  to  point  out  that  the  units  used  within  this  study 
are  often  mixed.  For  the  majority  of  the  study  the  SI  system  of  units  is  employed. 
However,  it  is  standard  in  the  literature  of  the  subject  to  use  units  of  feet  per  second  (fps) 
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for  velocity,  feet  (ft)  for  altitude,  and  nautical  miles  (nmi)  for  distance.  Where 
appropriate  these  units  will  be  given  in  addition  to  the  SI  units. 
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II,  Problem  Setup 


Chapter  Overview 

The  purpose  of  this  chapter  is  to  set  up  the  atmospheric  reentry  problem  for  the 
Space  Maneuver  Vehicle  (SMV).  First,  the  equations  of  motion  for  a  point  mass 
reentering  the  Earth’s  atmosphere  will  be  stated.  Next,  a  model  of  the  Space  Maneuver 
Vehicle  will  be  constructed.  Using  the  model  the  aerodynamic  properties  of  the  reentry 
vehicle  will  be  approximated  using  Newtonian  impact  theory.  Finally,  assumptions  will 
be  made  concerning  the  Earth’s  atmosphere  and  gravitational  field. 

Equations  of  Motion 

The  first  step  in  examining  a  body  reentering  the  Earth’s  atmosphere  is  to  develop 
a  set  of  equations  to  describe  the  motion.  The  set  of  equations  developed  and  used  herein 
is  similar  to  the  set  used  in  Vinh  as  well  as  in  Regan  and  Anandakrishnan.  The 
differences  are  that  the  flight  path  angle,  y,  will  be  positive  below  the  horizontal  plane, 
and  the  heading  angle,  \|/,  will  be  measured  from  North,  instead  of  East. 

The  coordinate  system  used  is  constructed  with  one  axis  aligned  with  the  Earth- 
relative  velocity  vector  (Harpold  and  Graves,  1978:106).  The  position  of  the  vehicle  is 
defined  by  a  vector  from  the  center  of  the  Earth  to  the  vehicle.  Then,  the  plane  formed  by 
this  vector  and  the  velocity  vector  is  perpendicular  to  the  second  axis,  with  the  third  axis 
completing  the  right-handed  coordinate  system  (Harpold  and  Graves,  1978:106).  The 
equations  or  motion  developed  in  this  coordinate  system  are: 
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V  =  -D  +  g  sin  Y 


(1) 


where; 


Vy  = 


f 

g - 


eosy-Zeoscj) 


V  eosyxj/  = 


V  y 


cos^  y  sin  \|/  tan  0  -  Z  sin  (j) 


V  =  Earth-relative  veloeity  (km/s) 

D  =  Drag  (km/s^) 

g  =  Gravitational  aeceleration  (km/s  ) 
y  =  Flight  path  angle 
r  =  (Re  +  h)  (km) 

Re  =  Radius  of  Earth  (km) 
h  =  Altitude  of  vehiele  (km) 

E  =  Eift  (km/s^) 

(j)  =  Bank  angle 
\|/  =  Heading  angle 
0  =  Eatitude 


(2) 

(3) 


The  derivation  of  these  equations  is  presented  in  Hypersonic  and  Planetary  Entry  Flight 
Mechanics  and  is  not  duplieated  herein  (Vinh,  1980:21-27).  It  is  worth  noting,  however, 
that  eertain  assumptions  have  been  employed  to  simplify  the  equations  to  the  form 
presented  above.  Assuming  a  non-rotating  Earth,  and  therefore  a  non-rotating 
atmosphere,  eliminates  the  co  and  co  terms  whieh  result  due  to  the  Earth’s  rotation  (Vinh, 
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1980:27).  In  other  words,  the  Coriolis  and  centripetal  acceleration  have  been  neglected. 
This  assumption  is  appropriate  since  the  effects  of  the  rotating  atmosphere  on  the  vehicle 
are  small  compared  with  aerodynamic  forces  due  to  the  vehicle’s  velocity  (Vinh,  1980:3). 
This  approach  is  consistent  with  the  methods  used  in  Shuttle  Entry  Guidance  (Harpold  & 
Graves,  1978:106). 

Figure  3  shows  the  coordinate  system  employed  in  this  analysis.  For  the  purposes 


Z 


Figure  3:  Reentry  Coordinate  System 

of  this  study  the  vehicle  is  considered  to  be  a  point  mass.  The  bank  angle,  (j),  is  defined  as 
positive  in  the  right-handed  sense  about  the  drag  vector  (ie.  opposite  the  velocity  vector). 
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That  is,  referencing  Figure  3,  positive  roll  is  defined  as  roll  to  the  left,  and  negative  roll  as 
roll  to  the  right.  The  heading  angle,  \|/,  is  measured  clockwise  from  North  to  the  velocity 
vector. 

The  set  of  equations  presented  (Eqns.  1-3)  are  the  kinetic  equations  of  motion  for 
this  system.  To  completely  describe  the  system,  the  kinematic  equations  must  also  be 
understood.  Once  again,  the  equations  used  in  this  study  are  similar  to  those  derived  in 
Hypersonic  and  Planetary  Flight  Mechanics  (Vinh,  1980:26).  The  kinematic  equations 
are: 


h  =  -V  sin  Y 
A  F  cosy  sin \|/ 

y  = - 

rcoscj) 

V  cos  Y  cos  \|/ 

k  =  - 

r 

where  A,  =  Longitude. 


(4) 

(5) 

(6) 


Vehicle  Model 

The  equations  developed  in  the  preceding  section  act  to  model  the  dynamics  of 
atmospheric  reentry.  The  purpose  of  this  study  is  to  examine  the  dynamics  of  a  particular 
lifting-body  reentry  vehicle — the  Space  Maneuver  Vehicle.  At  present,  the  SMV  exists 
in  concept  only,  however,  the  X-40A  and  X-37  programs  are  actively  developing 
technologies  in  support  of  its  development  (Shaw,  2000:1).  What  does  exist  for  the  SMV 
is  nominal  values  for  vehicle  parameters  ranging  from  weight,  to  length,  to  payload  bay 
dimensions  (Shaw,  2000:2). 
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In  order  to  study  the  reentry  characteristics  of  a  given  vehicle,  an  aerodynamic 
model  of  the  vehicle  must  be  developed.  The  approach  taken  in  this  study  is  to  build  a 
model  of  the  Boeing  SMV  concept  in  MATLAB  and  use  the  model  to  determine  the 
aerodynamic  characteristics  using  Newtonian  impact  theory.  Table  1  below  presents  the 
dimensions  of  the  SMV  taken  from  the  Boeing  SMV  concept: 


Table  1:  SMV  Parameters  (Shaw,  2000:2) 


Parameter 

Value 

Weight 

10000  lb 

Length 

29  ft 

Wingspan 

15  ft 

Height 

9.5  ft 

In  addition  to  the  parameters  presented  in  Table  1,  dimensions  including  nose  width. 


body  width,  and  chord  length  were  determined  based  on  interpretation  of  Boeing  X-37 


concept  drawings  (Figure  4). 


Figure  4:  Boeing  X-37  (X-37  2004) 
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These  values  were  then  used  to  ereate  a  60-panel  model  of  the  SMV  using  MATLAB,  see 
Figure  5.  A  more  detailed  diagram  is  shown  in  Appendix  A. 


Figure  5:  MATLAB  model  of  Spaee  Maneuver  Vehiele 


The  model  was  then  used  to  determine  the  forees  aeting  on  the  vehiele  through  the 
employment  of  Newtonian  impaet  methods.  Newtonian  impaet  theory  assumes  that  a 
partiele  of  fluid  whieh  impinges  upon  a  body  will  lose  all  momentum  normal  to  the 
surfaee,  while  the  tangential  momentum  will  be  unaffeeted  (Regan,  1993:352).  This 
method  ean  be  employed  to  determine  the  forees  aeting  on  a  flat  plate  (or  in  the  ease  of 
the  vehiele  model,  a  series  of  flat  plates)  in  both  the  free-moleeular  flow  and  hypersonie 
eontinuum  regions  (Regan,  1993:356).  The  FORT  AN  eode  used  for  this  analysis  was 
provided  by  Dr.  William  Wiesel,  Professor  of  Astronautieal  Engineering,  AFIT,  and  is 
similar  to  the  eode  provided  in  Dynamics  of  Atmospheric  Reentry  (Regan,  1993:321- 
328). 
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Aerodynamic  Properties, 

Upon  returning  from  orbit  the  vehiele  will  pass  from  the  near  vaeuum  eonditions 
of  spaee,  through  the  upper  regions  of  Earth’s  atmosphere,  to  the  denser  regions  of  the 
lower  atmosphere.  Throughout  the  reentry  trajeetory  the  vehiele  will  be  operating  in  one 
of  three  flight  regimes:  free  moleeular  flow,  transition,  or  eontinuum  flow  (Regan, 
1993:316).  The  regime  that  the  vehiele  is  operating  in  is  dependent  upon  the  Knudsen 
number.  The  Knudsen  number  is  a  non-dimensional  parameter  whieh  indieates  the 
relative  importanee  of  the  partieulate  nature  of  air  (Regan,  1993:3 16).  The  Knudsen 
number  ean  be  defined  as  (Regan,  1993:316): 

A,  mean  free  path 

Kn  =  —  = - - - - -  (8) 

L  characteristic  length 

The  mean  free  path  of  a  partiele  is  a  measure  of  the  relaxation  distanee  in  a  gas  (Regan, 
1993:316).  The  relaxation  distanee  is  a  measure  of  the  distanee  a  partiele  of  gas  will 
move  through  a  flowfield  before  interaeting  with  another  partiele  of  gas  (Regan, 

1993:3 13).  If  the  flowfield  has  a  high  number  density  of  partieles  then  the  eollision 
frequeney  will  be  high  and  this  distanee  will  be  short  (Regan,  1993:313).  In  this  ease,  the 
flowfield  ean  be  treated  as  a  continuum.  In  contrast,  in  a  region  where  the  number 
density  of  particles  is  low,  the  collisions  will  occur  less  frequently  and  the  relaxation 
distance  will  be  longer  (Regan,  1993:314).  In  this  case  the  flowfield  can  be  treated  as 
free  molecular  flow.  The  characteristic  length  is  typically  chosen  as  the  mean 
aerodynamic  chord  (MAC)  for  a  lifting  body. 

It  now  remains  to  define  what  is  meant  by  short  and  long  relaxation  distances. 

The  free  molecular  flow  region  is  defined  as  the  region  where  Kn  »  1 ,  while  the 
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continuum  region  is  defined  where  Kn  «  1  (Regan,  1993:316).  For  the  purposes  of  this 
study,  regions  where  Kn  >10  will  be  eonsidered  free  moleeular  flow  regions.  For  the 
region  where  Kn  <  .01 ,  the  flowfield  will  be  treated  as  a  eontinuum.  The  region  between 
the  two,  .01  <  Kn  <10  will  be  eonsidered  the  rarefied- flow  transition  region  (Blanehard, 
1994:550).  In  this  region  the  flow  eannot  be  treated  as  a  continuum,  nor  can  the 
partieulate  nature  of  the  gas  moleeules  be  negleeted.  In  this  region,  the  following 
empirieal  bridging  formula  will  be  used  (Blanehard,  1994:553): 

C^v  =exp[-0.2998 1(1 .3849- logio^Sin)'^''']  (9) 

//'logjoifn<  1.3849 
otherwise  Cn  =1.0 

C.4  =exp[-0.2262(l. 2042- logio  Ain)' (10) 
if  logjo  Kn<\  .2042 
otherwise  =1.0 

where  Cn  and  Ca  are  the  bridging  eoeffieients,  with  values  between  0  and  1.  These 
eoeflieients  are  then  used  to  bridge  the  transition  region  as  follows: 

Q=Q^+(Q^-QJQ  (11) 

C,=C^+(C,^-C,^)^  (12) 

where: 

Cn  =  Normal  aerodynamie  foree  eoeffieient 
Ca  =  Axial  aerodynamie  foree  eoeffieient 
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and  the  subscripts  “F”  and  “C”  denote  the  free  molecular  flow  and  continuum  regions, 
respectively. 

For  a  given  reentry  vehicle,  and  therefore  a  given  characteristic  length,  the 
Knudsen  number  as  a  function  of  altitude  can  be  determined.  The  Knudsen  number  can 
be  represented  as  (Regan,  1993:315): 


Kn  = 


m  /  {-4^71X5^ PqC  '''^) 
L 


(13) 


where: 

m  =  particle  mass  (kg) 
a  =  effective  diameter  of  gas  particles  (km) 
po  =  earth  surface  density  (kg/km  ) 

H  =  scale  height  (km) 

Figure  6  is  plot  of  Knudsen  Number  as  a  function  of  altitude,  with  the  mean  aerodynamic 
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Figure  6:  Altitude  vs.  Knudsen  number,  MAC=5.85ft 
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chord  (MAC=5.85  ft)  calculated  from  the  MATLAB  model  of  the  SMV.  Therefore,  the 
transition  region  from  free  moleeular  flow  to  hypersonie  continuum  flow  is  from 
approximately  130  km  to  70  km  altitude.  For  the  purposes  of  this  study,  reentry  begins  at 
400,000  ft  altitude,  or  approximately  122.2  km.  Therefore,  the  vehicle  will  begin  reentry 
in  the  transition  region  and  enter  the  continuum  region  at  approximately  70  km. 

To  this  point,  a  flat-panel  model  of  the  vehicle  has  been  created  and  a  method  for 
accounting  for  the  partieulate  nature  of  the  flow-field  has  been  developed.  The  goal  is  to 
determine  the  aerodynamic  coefficients  of  the  reentry  vehicle  for  use  in  the  equations  of 
motion  and  the  guidance  algorithm.  In  particular,  a  polynomial  function  of  angle  of 
attaek,  a,  is  desired.  From  the  aforementioned  Newtonian  impact  theory  analysis  the  lift 
and  drag  coefficients  as  functions  of  angle  of  attack  for  both  the  free  molecular  flow 
region  and  the  hypersonic  continuum  flow  region  have  been  developed.  The  resulting  lift 
coefficients  are  shown  in  Figures  7  and  8. 
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Figure  7:  Hypersonic  Continuum  Lift  Coefficient 
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Figure  8  shows  the  expeeted  result  that  very  little  lift  is  being  generated  in  the  free 
moleeular  flow  region. 


AOA  in  degrees 

Figure  8:  Free  Moleeular  Flow  Lift  Coeffieient 
Figures  9  and  10  show  the  determined  drag  eoeffieients  in  the  hypersonie  eontinuum  and 
free  moleeular  flow  regions,  respeetively. 
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Figure  9:  Hypersonie  Continuum  Drag  Coeffieient 
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Figure  10:  Free  Molecular  Flow  Drag  Coefficient 


To  determine  polynomial  expressions  for  the  aerodynamic  coefficients  in  each 
region,  a  least  squares  method  was  employed.  First,  it  was  determined  that  a  quadratic 
expression  of  the  aerodynamic  coefficient  in  terms  of  angle  of  attack  would  be  used. 

Next,  since  reentry  at  a  high  angle  of  attack  is  desirable  to  minimize  aerodynamic 
heating,  the  range  of  angle  of  attack  over  which  the  vehicle  would  operate  during  reentry 
was  established  as  15°-50°  (Chapman,  1967:37).  For  each  case  (ie.  lift  coefficient  in  the 
continuum  region),  this  defined  an  over-determined  set  of  equations  in  three  unknowns. 
This  system  of  equations  can  be  expressed  as: 

[A\x  =  bLc  (14) 

The  matrix  [^4]  has  three  columns.  The  first  column  is  the  square  of  the  angle  of  attack, 
the  second  column  is  the  angle  of  attack,  and  the  third  column  is  1 .  Thus: 
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(15) 


AOA^^ 

AOA^ 

AOAy 

AOA„ 

The  elements  of  the  veetor  bic  are  the  lift  eoefficients  in  the  hypersonic  continuum 
region  corresponding  to  the  particular  angle  of  attack.  Thus: 


^LC  ~ 


C 


LC\ 


c 


LCn 


(16) 


Since  the  system  of  equations  is  overdetermined,  it  is  unlikely  that  there  will  be  a 
solution,  X ,  which  exactly  fits  the  data  in  .  The  least  squares  solution  to  an 
overdetermined  system  such  as  this  satisfies: 

lAnA]X  =  [Af  bic  (17) 


where  the  vector,  X ,  contains  the  coefficients  of  the  quadratic  (Strang,  1988:156): 

a(AOA.y+b(AOA.)  +  c  =  C,c,  (18) 

Since  the  columns  of  the  matrix  [^4]  are  linearly  independent,  the  matrix  [^4]^  [^4]  is 

invertible,  and  the  solution  to  the  least  squares  problem  can  be  written  as  (Strang, 
1988:156): 

X  =  {[A\^[Ay[A\bLc  (19) 

Then,  solving  for  the  coefficients,  the  lift  coefficient  in  the  hypersonic  continuum  flight 
regime  is  approximated  by: 

Qc  =-3-3069*10  "(^0^)"+0.0552(^0^)-0.5034  (20) 
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The  preceding  procedure  is  repeated  to  determine  the  lift  coefficient  in  the  free-molecular 


flow  region  simply  by  replacing  the  vector,  bic ,  with  the  vector,  biF ,  which  is  made  up 
of  the  lift  coefficients  in  the  free  molecular  flow  region.  Then  the  lift  coefficient  in  the 
free  molecular  flow  region  is  approximated  by; 

=-3.8241*10  '(^0^)"+0.0033(^0^)-0.0022  (21) 

Figures  1 1  and  12  show  the  polynomial  approximation  of  the  lift  coefficients  over  the 
region  of  interest  in  the  hypersonic  continuum  region  and  the  free  molecular  flow  region, 
respectively. 


Figure  11;  Quadratic  Expression  of  Hypersonic  Continuum  Lift  Coefficient 

The  drag  coefficients  in  each  region  of  flight  are  determined  by  the  same  procedure 
outlined  above.  The  resulting  polynomial  approximations  of  the  drag  coefficients  are; 

=-9-6602*10^(HOH)'-0.0150(HOH)  +  0.2402  (22) 
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Figure  12;  Quadratic  Expression  of  Free  Molecular  Flow  Fift  Coefficient 


for  the  hypersonic  continuum  region,  and 

=  -1 .3241  *  10  +  0.0373(^0^)  +  0.2265 

for  the  free  molecular  flow  region.  Figures  13  and  14  show  the  polynomial 


Figure  13;  Quadratic  Expression  of  Flypersonic  Continuum  Drag  Coefficient 


(23) 
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Figure  14;  Quadratic  Expression  of  Free  Molecular  Flow  Drag  Coefficient 
approximations  to  the  drag  coefficients  in  the  hypersonic  continuum  and  free  molecular 
regions,  respectively. 

Physical  Assumptions 

For  this  analysis  the  gravitational  field  of  the  earth  is  assumed  to  be  constant. 

That  is,  the  magnitude  of  the  gravitational  force  is  independent  of  altitude  and  position. 
The  atmosphere  is  modeled  as  an  exponential  atmosphere  described  by  (Regan,  1993:38): 

p  =  Po  exp(-/z  /  H)  (24) 

where; 

Po  =  sea-level  density  (1.752*10^9  kg/km^3) 

H  =  Scale  height  (6.7  km) 

This  atmospheric  model  shows  excellent  agreement  with  the  1976  Standard  Atmosphere 
to  an  altitude  of  approximately  120  km  (Regan,  1993:39). 
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Ill,  Methodology 


Chapter  Overview 

The  purpose  of  this  ehapter  is  to  develop  the  guidanee  approaeh  that  will  be  used 
to  guide  the  SMV  during  reentry.  First,  a  set  of  initial  and  final  eonditions  whieh  bound 
the  trajeetory  will  be  developed.  Next,  the  eoneept  of  drag  profile  guidanee  will  be 
developed  and  applied  to  this  partieular  vehiele.  Finally,  the  guidanee  and  eontrol 
algorithms  will  be  developed. 

Boundary  Values 

The  guidanee  approaeh  used  will  guide  the  reentry  vehiele  from  the  eonditions  at 
an  initial  state  to  the  eonditions  at  a  final  state  without  violating  any  reentry  eonstraints. 
In  this  seetion,  the  initial  and  final  eonditions  will  be  developed. 

In  a  somewhat  eounter-intuitive  approaeh,  the  final  eonditions  will  first  be 
developed,  and  then  the  initial  eonditions.  This  approaeh  is  driven  by  the  faet  that  a  key 
element  of  a  sueeessful  reentry  is  that  the  vehiele  reaehes  the  desired  termination  point  at 
an  aeeeptable  energy  state.  Therefore,  the  final  position  is  fixed  to  ensure  this  eondition 
is  met,  and  then  the  initial  eonditions  are  determined  based  upon  range  eonstraints. 

The  reentry  guidanee  for  the  Spaee  Shuttle  is  designed  to  “eontrol  the  entry 
trajeetory  from  initial  penetration  of  the  Earth’s  atmosphere  (altitude  of  400,000  ft  and 
range  of  approximately  4100  nmi  from  runway)  until  aetivation  of  the  terminal  area 
guidanee  at  an  Earth-relative  speed  of  2500  fps”  (Harpold  and  Graves,  1978:99).  This 
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paper  will  not  explore  the  guidanee  that  controls  the  vehicle  to  the  runway  from  the 
termination  point. 

To  determine  the  location  of  the  termination  point,  the  value  of  Earth-relative 
velocity  of  2500  fps  (0.764  km/s)  used  for  the  Space  Shuttle  will  be  used  as  the 
terminating  velocity.  The  altitude  at  termination  of  entry  guidance  will  be  29.4  km 
(96,250  ft)  and  the  distance  from  the  termination  point  to  the  runway  will  be  30  km  (16.1 
nmi).  The  range  that  the  vehicle  is  capable  of  flying  from  this  point  can  be  estimated  by 
(Harpold  and  Graves,  1978:109); 


R^iE,  Ec)in^ 


(25) 


where; 

2  2 

Ef  =  specific  energy  at  termination  point  (km  /s  ) 

2  2 

Eg  =  specific  energy  at  landing  (km  /s  ) 

Dp  =  drag  at  termination  point  (km/s  ) 

Dg=  drag  at  landing  (km/s  ) 

With  the  values  above  for  termination  velocity  and  altitude,  and  a  landing  speed 
of  250  fps,  the  remaining  range  of  the  vehicle  is  approximately  64  km  (34.5  nmi).  This 
excess  energy  would  either  be  dumped  through  wide  turns  or  used  to  align  the  vehicle 
with  the  runway.  An  alternate  method  of  predicting  the  range  to  go  from  the  termination 
point  based  upon  integration  of  the  equations  of  motion  is  presented  in  Appendix  B.  The 
two  methods  show  agreement,  with  the  method  presented  herein  being  the  more 
conservative. 
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By  defining  the  entry  guidanee  termination  values  for  veloeity  and  altitude,  the 
terminal  value  for  drag  is  also  defined,  since; 

D  =  ^pV^C^S  (26) 

and  the  angle  of  attack  at  termination,  and  therefore  the  drag  coefficient,  is  specified. 

Next,  the  physical  location  on  the  Earth  of  the  entry  guidance  termination  point 
will  be  determined.  For  this  paper,  the  landing  site  is  assumed  to  be  at  Kennedy  Space 
Center  (KSC),  FL.  The  coordinates  of  KSC  are; 

Fatitude;  28.6°  N 
Fongitude;  279.4°  E 

The  entry  guidance  termination  point  will  be  determined  to  lie  30  km  directly  West  of 
KSC.  Since  one  degree  of  longitude  (°F)  is: 

°F=2;rRg  cos{Latitude)  /  360 
the  coordinates  of  the  entry  guidance  termination  point  are; 

Fatitude;  28.6°  N 
Fongitude;  279.1°  E 

Table  2  summarizes  the  conditions  at  entry  guidance  termination; 


Table  2:  Terminal  Conditions  for  Entry  Guidance 


Conditions 

Value 

Velocity 

0.764  km/s 

Altitude 

29.4  km 

Fatitude 

28.6°  N 

Fongitude 

279.1°  E 
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Next,  the  initial  conditions  for  the  reentry  trajectory  must  be  determined.  The 
reentry  trajectory  will  begin  at  an  altitude  of  122.2  km  (400,000  ft) — the  same  altitude 
used  for  Shuttle  entry  guidance  (Harpold  and  Graves,  1978:99).  The  velocity  at  reentry 
will  be  dependent  upon  the  altitude  of  the  vehicle  orbit.  For  this  study,  the  initial  orbit  of 
the  vehicle  will  be  assumed  circular.  The  de-orbit  burn  will  be  the  first  bum  of  a 
Hohmann  transfer  from  the  initial  orbit  to  a  circular  orbit  of  80  km,  as  illustrated  in 
Figure  15. 


Figure  15:  De-Orbit  and  Reentry  Interface 
For  the  transfer  ellipse  (Wiesel,  1997:56): 

a^  +  af 

a,  = - - 

'  2 

and 

rp=at(\-e) 

where  rp  is  the  perigee  point  of  the  transfer  ellipse,  e  is  the  eccentricity  and: 
at  =  semi-major  axis  of  transfer  ellipse  (km) 


(27) 


(28) 
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ai  =  semi-major  axis  (radius)  of  initial  orbit  (km) 
af  =  semi-major  axis  (radius)  of  80  km  destination  orbit 
The  equations  above  can  be  solved  to  find  the  semi-major  axis  and  eccentricity  of  the 
transfer  ellipse.  These  quantities  can  then  be  used  to  determine  the  magnitude  of  the 
angular  momentum,  h,  from  (Wiesel,  1997:56): 

h  =  ^J\Jia(l-e^)  (29) 

Then,  letting  r  be  equal  to  the  distance  from  the  center  of  the  Earth  to  the  vehicle  at  the 
reentry  altitude  (ie.  r  =  Re  +  122.2  km),  the  initial  velocity  can  be  found  from  (Wiesel, 
1997:74): 


\  r  a 


(30) 


Then  the  initial  flight  path  angle,  Yq  ,  can  be  determined  from  (South,  1974:159): 


Yo  =  cos 


-1  h 


rVr. 


(31) 


The  preceding  equations  have  specified  the  initial  energy  state  for  the  reentry 
vehicle.  Now  the  location  of  the  vehicle  at  the  beginning  of  the  reentry  trajectory,  with 
respect  to  a  non-rotating  earth,  must  be  determined.  The  position  of  the  vehicle  at  the 
initiation  of  the  reentry  trajectory  will  be  a  function  of  the  projected  range  of  the  vehicle. 
This  projected  range,  in  turn,  will  be  a  function  of  the  vehicle’s  initial  velocity,  pre¬ 
determined  drag  profile,  and  terminal  conditions.  The  drag  profile  which  the  vehicle  will 
follow  will  be  developed  in  the  next  section.  As  such,  further  discussion  of  the  position 
of  the  vehicle  at  reentry  is  presented  in  Appendix  C.  The  initial  conditions  for  reentry 
from  an  initial  circular  orbit  of  500  km  are  summarized  in  Table  3. 
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Table  3:  Initial  Conditions  for  Entry  Guidance  (ho=500km) 


Conditions 

Value 

Velocity 

7.940  km/s 

Altitude 

122.2  km 

Fatitude 

Orbit  Dependent 

Fongitude 

Orbit  Dependent 

Drag  Profile  Guidance 

This  section  will  develop  the  concept  of  controlling  the  vehicle’s  reentry 
trajectory  by  controlling  to  a  desired  drag  profile.  First  the  concept  of  the  drag-velocity 
(D-V)  plane  and  the  limiting  constraints  will  be  developed.  Next,  a  three-phase  nominal 
reentry  profde  will  be  developed  in  the  drag-velocity  plane. 

Reentry  Constraints. 

For  a  vehicle  reentering  the  Earth’s  atmosphere  there  are  four  primary  constraints 
that  determine  the  reentry  corridor.  For  the  purposes  of  this  study,  these  constraints  will 
be  plotted  in  the  D-V  plane.  During  the  early  phase  of  entry  the  vehicle  will  be  subjected 
to  high  heating  rates  as  it  descends  into  the  increasingly  dense  atmosphere.  While 
heating  rates  are  the  primary  concern  during  the  early  phase  of  reentry,  later  in  flight 
vehicle  normal  loads  and  dynamic  pressure  loads  become  increasingly  important.  The 
final  constraint  is  the  range  constraint.  A  functional  guidance  program  will  guide  a 
vehicle  through  the  reentry  corridor  without  violating  the  physical  constraints  while 
delivering  the  vehicle  to  a  pre-determined  termination  point  at  an  acceptable  energy  state. 
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Heating  Rate, 

For  a  vehicle  reentering  the  earth’s  atmosphere  two  heat  related  quantities  must  be 
accounted  for:  heating  rate  and  heating  load.  For  a  vehicle  that  uses  an  ablation  type 
heat  shield,  the  total  heat  load,  or  heat  absorbed,  is  of  primary  importance  (Chapman, 
1967:38).  For  a  radiation-cooled  vehicle  the  maximum  heating  rate  is  the  driving  factor 
(Chapman,  1967:38).  As  mentioned  earlier,  the  SMV  exists  only  in  concept  at  this  point, 
however,  the  Boeing  X-37  will  test  possible  thermal  protection  system  (TPS)  materials 
which  could  be  used  on  an  operational  SMV  (Shaw  2000:1 14).  Therefore,  for  the 
purposes  of  this  study,  the  assumption  will  be  made  that  the  SMV  reentry  heating 
constraint  will  be  controlled  by  the  heating  rate. 

The  heating  rate  constraint  equation  is  (Lu,  1997:145): 

ry~  4^0  V  ^  gmax  (32) 

\  Po.re/  _  Ve/  _  _  4^eSo  _ 

where: 

po,ref=  1  kg/m^ 

Vref=  1  rn/s 

^max  maximum  stagnation  point  heating  rate  (W/m  ) 

=  heat  flux  transmission  coefficient  (W/m  ) 

The  parameters  for  maximum  stagnation  point  heating  rate  ( )  and  heat  flux 
transmission  coefficient  ( )  are  taken  from  Lu  for  a  reusable  launch  vehicle  (RLV)  and 

based  on  a  reference  sphere  of  Im  (Lu,  1997:145): 

^.ax=  544,300  W/m' 
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2 


C^=  1.65*10-^  WW 


Then,  Equation  (32)  can  be  rearranged  to  yield; 


max 


(33) 


This  expression  can  then  be  substituted  into  the  drag  equation  to  give; 

D  1  PoEe/  ^max  ^ 

-  =  ^C^S  \  (34 

m  2m  V 

All  that  remains  to  do  to  quantify  the  heating  constraint  is  to  specify  the  angle  of 
attack  profile.  The  angle  of  attack  profile  is  chosen  such  that  the  reentry  vehicle  initially 
enters  the  atmosphere  with  a  high  angle  of  attack  and  maintains  that  angle  of  attack  until 
a  specified  velocity.  At  the  specified  velocity,  the  angle  of  attack  will  decrease  linearly 
to  the  angle  of  attack  corresponding  to  max  L/D  at  the  final  conditions.  Figure  16  shows 
the  angle  of  attack  profile  used  in  this  study  for  the  SMV. 
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Figure  16;  Angle  of  Attack  profile 
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Using  the  angle  of  attaek  profile,  the  drag  per  unit  mass  ean  be  plotted  versus  veloeity  as 
shown  in  Figure  17. 
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Figure  17:  Ideating  Constraint  in  D-V  Plane 


It  should  be  noted  that  drag  is  plotted  against  negative  veloeity  as  a  eonvention. 

In  the  D-V  plane  the  vehiele  will  be  traveling  from  left  to  right  as  it  enters  the  atmosphere 
and  deeelerates  to  its  terminal  point.  The  heating  eonstraint  represents  the  first  eonstraint 
to  reentry — developing  the  other  eonstraints  will  develop  the  reentry  eorridor. 


Maximum  Normal  Load. 

The  maximum  normal  load  eonstraint  represents  the  maximum  load  eonstraint  in 
the  body-normal  direetion  (Lu,  1997:145).  Figure  18  shows  the  relationship  of  the  forees 
aeting  on  a  reentry  body: 
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In  the  above  figure,  z  represents  the  body  normal  axis,  and  a  represents  the  longitudinal 
axis  of  the  vehicle.  As  can  be  seen  from  the  figure,  both  lift  and  drag  contribute  to  the 
normal  loads  on  the  vehicle.  From  the  geometry  it  can  be  seen  that  the  normal  load  on 
the  vehicle  is: 


Substitute, 


then. 


Nz  =— — cosaH — sma 
m  m 


(35) 


(36) 


(37) 


Let  the  maximum  normal  load  constraint  be  Az.max .  Then,  rearrange  and  solve  for  D/m: 


D  ^  Az,max 

m  . 

— ^cosa  +  sma 


(38) 
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From  Harpold  and  Graves,  the  maximum  normal  load  boundary  for  the  Spaee  Shuttle  is 


defined  to  be  2.5g  (1978:103).  Therefore,  setting  A^z.max  =  2.5g  and  plotting  in  the  D-V 
plane  yields  the  normal  load  eonstraint  as  shown  in  Figure  19. 


Dynamic  Pressure  Constraint, 

The  dynamic  pressure  constraint  to  reentry  can  be  seen  directly  from  the  drag 
equation.  Letting  qmaxbe  the  maximum  allowable  dynamic  pressure  leads  directly  to: 


D  C^S 

—  ^^max  - 

m  m 


(39) 


Using  the  Space  Shuttle  value  of  qmax  =  16,280  N/m^  and  plotting  in  the  D-V  plane  yields 
the  dynamic  pressure  constraint,  as  shown  in  Figure  20  (Flarpold  and  Graves,  1978:103): 


Figure  19:  Normal  Load  Constraint  in  D-V  Plane 
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Figure  20;  Dynamic  Pressure  Constraint  in  D-V  Plane 


Equilibrium  Glide  Constraint, 

The  final  constraint  to  determine  in  the  D-V  plane  is  the  equilibrium  glide 
constraint.  The  equilibrium  glide  constraint  is  associated  with  the  maximum  range 
achievable.  To  determine  the  constraint,  the  flight  path  angle  is  held  constant  and  the 
bank  angle  is  set  to  0,  therefore  achieving  the  minimum  drag  level  (Lu,  1997:145).  Then, 
from  Equation  (2)  with  the  flight  path  angle  assumed  small,  ie  y  =  0 ; 


-H— 

m  Re-l- A 


+  g) 


(40) 


Substitute; 


and. 


m  m  D 


(41) 


35 


(42) 


D 


to  get; 


D  C  -V^ 
m  C,  Re+  h 


(43) 


Plotting  this  equation  in  the  D-V  yields  the  equilibrium  glide  slope  constraint  as  shown  in 
Figure  21; 
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Figure  21;  Equilibrium  Glide  Slope  Constraint 


Plotting  all  four  constraints  together  defines  the  reentry  corridor  in  the  D-V  plane 
as  shown  in  Figure  22.  An  effective  entry  guidance  controller  will  guide  the  vehicle’s 
trajectory  in  such  a  way  as  to  not  violate  any  of  the  determined  entry  constraints,  while  at 
the  same  time  controlling  the  vehicles  energy  such  that  the  entry  guidance  termination 
point  is  reached  at  the  desired  energy  state. 
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Figure  22;  D-V  Plane  Reentry  Corridor  for  SMV 


Reference  Profile 

The  previous  section  defined  the  reentry  corridor  in  the  D-V  plane  for  the  SMV. 
Next  a  reference  profile  that  falls  within  the  corridor  will  be  developed.  The  reference 
profile  will  be  constructed  of  three  segments;  the  temperature  control  segment,  the 
constant  drag  segment,  and  the  transition  segment  (Flarpold  &  Graves,  1978;  109). 

Temperature  Control , 

The  primary  concern  during  the  early  stages  of  atmospheric  reentry  is  the  heating 
rate  constraint.  The  reference  profile  is  developed  by  choosing  a  drag  acceleration  level 
at  the  initial  velocity  which  lies  within  the  reentry  corridor.  Next,  a  termination  velocity 
is  chosen  for  the  temperature  control  segment.  For  the  termination  velocity,  a  drag 
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acceleration  level  within  the  reentry  eorridor  is  chosen.  Then,  drag  acceleration  at  any 
point  in  the  temperature  eontrol  region  is  given  by: 

D  =  q+C2V  +  qV^  (44) 

where  the  eonstant  C3  =  3*10'"^  km'^  and  Ci  and  C2  are  determined  by  the  initial  and  final 
values  of  velocity  and  drag  (Harpold  &  Graves,  1978:109). 

Constant  Drag, 

The  initial  velocity  and  drag  acceleration  for  the  eonstant  drag  segment  are  the 
same  as  the  termination  velocity  and  drag  acceleration  for  the  temperature  control  region. 
Sinee  the  drag  aeeeleration  level  remains  constant  throughout  this  region, 

B  =  C,  (45) 

defines  the  drag  level  (Harpold  &  Graves,  1978:109).  Then,  a  termination  velocity 
which,  given  the  drag  level,  lies  within  the  reentry  eorridor  is  ehosen. 

Transition, 

The  transition  region  controls  the  vehiele  from  its  constant  drag  profile  to  the 
desired  conditions  at  entry  termination.  The  initial  eonditions  for  the  transition  region 
will  be  the  same  as  the  final  eonditions  of  the  constant  drag  region.  The  termination 
conditions  are  dictated  by  the  boundary  conditions  for  entry  guidance  termination 
established  earlier.  Given  those  eonditions,  the  following  equation  ean  be  solved  for  C5 
and  the  drag  acceleration  level  at  any  point  in  the  region  can  be  determined  (Harpold  & 
Graves,  1978:109): 
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D  —  Dp  +  Cj  (£  —  Ep  ) 


(46) 


where; 


Dp  =  drag  aeeeleration  level  at  entry  guidanee  termination  (km/s^) 

2  2 

Ep  =  speeifie  energy  at  entry  guidanee  termination  (km  /s  ) 

The  referenee  profile  ean  then  be  plotted  in  the  D-V  plane  with  the  reentry  eorridor  as 
shown  in  Figure  23; 


Figure  23;  D-V  Plane  Referenee  Profile  for  SMV 


Guidance 

The  guidanee  approaeh  applied  here  is  the  approaeh  used  for  Shuttle  entry 
guidanee.  The  guidanee  algorithm  “is  designed  on  the  prineipal  of  analytieally  defining  a 
desired  drag  aeeeleration  profile  and  eommanding  the  (vehiele)  attitudes  to  aehieve  the 
desired  profile”  (Harpold  and  Graves,  1978;  103).  With  the  referenee  drag  profile 
defined,  it  remains  to  determine  a  method  to  eontrol  the  vehiele  to  this  profile.  The  two 
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parameters  available  for  eontrolling  the  vehicle  to  this  profile  are  angle  of  attack,  a,  and 
bank  angle,  (j).  Using  bank  angle  as  the  primary  trajectory  control  parameter  allows  the 
selection  of  an  angle  of  attack  profile  which  will  provide  desirable  geometry  for  thermal 
protection  (Harpold  &  Graves,  1978:105).  Thus,  the  bank  angle  is  modulated  to  control 
the  vehicle  to  the  reference  drag  profile.  In  doing  this,  the  total  entry  range  and 
crossrange  are  also  controlled  by  the  bank  angle  (Harpold  &  Graves,  1978:105).  Then, 
angle  of  attack  modulation  is  used  for  short  period  trajectory  control,  such  as  during  a 
roll-reversal  when  the  bank  angle  will  roll  through  0°  and  thus  drive  the  vehicle  to  a 
lower  than  desired  drag  acceleration  level. 

Range  Modulation, 

The  guidance  algorithm  is  implemented  at  1  sec  intervals  throughout  the  reentry 
trajectory.  At  each  instance,  the  reference  drag  acceleration  profile  is  updated  to  account 
for  differences  between  nominal  range  to  go  and  actual  range  to  go.  The  nominal  range 
to  go  is  a  summation  of  the  nominal  range  to  go  in  each  of  the  three  phases  of  reentry- 
temperature  control,  constant  drag,  and  transition — and  is  calculated  by  (Harpold  and 
Graves,  1978:109): 
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Q  =  4C,q-C,^ 
ifQ>0 


-1  C  +C  V  +CV 


2C 


3 


Ci+QF  +C3F' 


C  ,  i2CF., +C,  i,2CF+C,, 

+ - W[tan^  [ — — ^]-tan  ‘[ — ^  ^]] 


c,^/e 


4q 


^/e 


,ye<0 


-1  C+CF  +CF 
i?,  =  —  ln[^  ^  ^ 

‘  2C3  C1+C2F+C3F' 


C 

+  — 


2C,V,,+C,-^ 

2C^Vpi  +C2+ 

2C,V+C,-^ 

2C,V+C,+4^ 


(47) 


where; 


Ri  =  Range  to  go  in  the  temperature  eontrol  region  (km) 

V  =  eurrent  vehiele  veloeity  (km/s) 

Vfi  =  veloeity  at  termination  of  temperature  eontrol  region  (km/s) 

Nominal  range  ealeulation  in  the  eonstant  drag  region  is  done  using  (Harpold  and 
Graves,  1978:109): 


where: 


R,= 


-V  2 

'  'FI 

2Q 


(48) 


Vf2  =  veloeity  at  termination  of  eonstant  drag  region  (km/s) 
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For  the  transition  region,  the  nominal  range  to  go  is  ealeulated  by  (Harpold  and  Graves, 
1978:109); 


^3  = 


(D-D,,) 


In 


D 


D 


F2 


where: 

2  2 

Ef3  =  speeifie  energy  at  termination  of  entry  guidanee  (km  /s  ) 

2 

Df3  =  drag  aeeeleration  level  at  termination  of  entry  guidanee  (km/s  ) 


(49) 


During  the  first  phase  of  reentry,  only  the  temperature  eontrol  referenee  drag 
aeeeleration  profile  is  adjusted — the  eonstant  drag  and  transition  referenee  profiles  are 
not  altered  (Harpold  and  Graves,  1978;  1 19).  The  effeet  of  adjusting  the  referenee  profde 
in  this  manner  is  “to  drive  the  trajeetory  baek  to  nominal  at  the  start  of  the  eonstant  drag 
phase”  (Harpold  and  Graves,  1978;  1 19).  The  same  approaeh  is  taken  during  the  eonstant 
drag  phase  where  the  transition  referenee  drag  profile  remains  unehanged. 

Next,  the  aetual  range  from  the  vehiele  to  the  termination  point  is  determined. 

The  range,  or  rather  the  range  angle,  0,  ean  be  determined  using  the  eurrent  latitude  and 
longitude  of  the  vehiele  and  the  latitude  and  longitude  of  the  termination  point,  from 
(Pratt,  2003:34): 

oos0  =  oos(T^,)  oos(T^  )cos{l^-lj-)  +  sin(T^, )  sin(T^  )  (50) 

where; 

Lv  =  current  latitude  of  vehicle 
Lf  =  latitude  of  termination  point 
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ly  =  current  longitude  of  vehiele 
If  =  longitude  of  termination  point 


Then,  the  aetual  range  to  go  is  determined  by: 

R  =  Re0 

and  the  range  error,  AR, 

AR  =  R-Rj, 

where: 


(51) 


(52) 


Rn  =  Nominal  range  to  go  (km)  -determined  by  Equations  (47-49) 

In  the  temperature  control  phase  of  reentry,  the  constant  Cl  will  be  determined  such  that 
the  updated  drag  reference  profile  meets  range  eonstraints.  Then  using  the  approximation 
for  the  derivative  of  range  with  respeet  to  drag  in  the  temperature  control  region 
(Tragesser  2003): 


dR^  +  7?2 

dD  ^f\  \ 

^  2  ^ 


Since, 


JRj  _  dR^  dC^  _  dC^ 
dD  dC^  dD  dD 

then,  substituting  and  solving  for  dCi: 


dC, 


(^D^+D^ 

R,  1 


(53) 


(54) 


(55) 


Then,  for  small  time  intervals. 
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(56) 


Ai?j  =  dR^ 

and  R1  is  calculated  from  Equation  (47)  above.  Then,  if  AR  is  negative  (nominal  range  is 
greater  than  aetual  range)  the  vehiele  will  be  guided  toward  a  higher  referenee  drag 
profde  sinee  Ci  will  be  inereased  by  dCi. 

As  the  reentry  progresses  the  vehiele  will  transition  from  the  temperature  eontrol 
phase  to  the  eonstant  drag  phase.  Once  again,  only  the  drag  profile  within  this  phase  will 
be  adjusted  to  null  range  errors.  In  the  constant  drag  phase,  the  eonstant  C4  is  modified  to 
adjust  the  drag  referenee  profile.  The  variation  of  range  with  C4  is  immediately  apparent 
from; 


dR^  dR^  dC,  dR-, 

— ^  = ^ ^  = ^  (57) 

dD  dC^  dD  dC^ 

since  D=C4.  Then,  with  dR2  =  AR,  differentiating  Equation  (48)  and  solving  for  dC4 
gives  (Tragesser  2003); 


dC, 


-IC/dR^ 

-V  2 

'  'FI 


(58) 


As  in  the  temperature  eontrol  region  if  AR  is  negative,  dC4  will  be  positive,  eausing  the 
vehiele  to  be  guided  toward  a  higher  drag  profde. 

Eor  the  transition  phase  of  reentry,  the  eonstant  C5  will  be  modified  to  null  range 
errors.  The  derivative  of  range  with  respeet  to  C5  ean  be  determined  by  first  rearranging 
the  drag  equation  for  the  transition  phase  (Equation  46)  and  solving  for  C5; 


_(D-Df3) 

(E-E,,) 


and  substituting  into  Equation  (49); 


(59) 
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(60) 


R.  =  —  In 
C 

'-5 


D 


D 


F3 


Then,  taking  the  derivative  with  respeet  to  C5  yields  (Tragesser  2003): 


JC  =■ 


dR. 


a 


-R.  +  ^ 


D 


(61) 


Onee  again,  with  dR3  =  AR,  C5  will  be  inereased,  and  therefore  the  drag  profde  inereased 
when  the  nominal  range  exeeeds  the  aetual  range. 


Bank  Angle  Modulation. 

The  previous  seetion  determined  the  varianee  of  Ci,  C4,  and  C5  to  aeeount  for 
range  errors.  This  seetion  will  show  how  these  updated  eoeffieients  are  used  to 
determine  the  appropriate  angle  of  attaek  to  guide  the  vehiele  to  the  referenee  drag 
aeeeleration  profile. 

Onee  the  updated  drag  aeeeleration  profde  is  defined,  the  desired  L/D  and  altitude 
rate  eorresponding  to  the  drag  aeeeleration  profile  need  to  be  determined  (Harpold  and 
Graves,  1978:109).  These  parameters  are  then  used  in  the  eontrol  law  to  eommand  the 
vehiele  to  the  referenee  profile  (Harpold  and  Graves,  1978:109).  The  equations  for  L/D 
and  altitude  rate  are  presented  in  Shuttle  Entry  Guidance  and  listed  in  Table  4,  with  the 
approximation  thatC^  =  =  0  (Harpold  and  Graves,  1978:109).  Using  these 

equations,  a  “eontrol  law  based  on  linearized  analysis  of  the  flight  dynamies,  whieh 
ensures  damping  of  oseillatory  type  trajeetory  motion”  ean  be  implemented  (Harpold  and 
Graves,  1978:113). 
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Table  4:  Reference  Trajectory  Parameters 


The  resulting  control  law,  with  f4=0  is  (Harpold  and  Graves,  1978:1 13): 

(^l=(^).+^(0-0.)  +  /,(*-Ao)  (62) 

where: 

(— ),  =  commanded  L/D 
D 

The  theoretical  gains  (fi  and  fi)  in  the  control  law  are  complex  and  lengthy  equations,  and 
therefore,  are  not  used  in  the  flight  software  (Harpold  and  Graves,  1978:123).  Instead, 
the  following  empirical  curve  fits  as  functions  of  drag  are  used  to  determine  the  gains: 
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(63) 


f,=fP^^+f,{D-D,) 

h  = 

where  the  coefficients  “fs,  fe,  f?,  fg,  and  fg  are  empirically  derived  constants  used  to  fit  the 
fi  and  f2  theoretical  functions”  (Harpold  and  Graves,  1978:123).  In  this  study,  the  curves 
used  for  f  and  f2  for  Shuttle  entry  guidance  were  used  to  determine  the  empirical 
coefficients.  In  the  case  of  the  controller  gain,  f2,  the  values  used  herein  are  the  negatives 
of  those  provided  in  Harpold  and  Graves — this  is  attributable  to  the  definition  of  positive 
flight  path  angle  below  the  horizontal  (1978:124).  Figures  24  and  25  show  the  controller 
gains  as  functions  of  drag. 
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Figure  24:  Controller  Gain  (f) 
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Figure  25;  Controller  Gain  (f2) 


The  (L/D)c  determined  from  the  eontrol  law  is  then  used  to  determine  the  bank 
angle  from  (Flarpold  and  Graves,  1978:124); 


(j)^  =  cos 


(L/Dl 

{LID) 


+  /i(a-ao) 


(64) 


where; 


(t)c  =  commanded  angle  of  attack 
L/D  =  current  lift-to-drag  ratio 
fi  1  =  controller  gain 
a  =  current  angle  of  attack 
ao  =  reference  angle  of  attack 
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In  addition,  angle  of  attack  modulation  is  used  for  short  period  drag  control,  such 
as  is  required  during  a  roll  reversal  (Harpold  and  Graves,  1978:125).  Angle  of  attack 
modulation  is  controlled  by  (Harpold  and  Graves,  1978:125): 

Aa  =  a-ao  =C^(Z)o-Z))/y;o  (65) 

where  fio  is  a  controller  gain.  During  a  roll  reversal,  the  angle  of  attack  will  be  increased, 
thus  creating  higher  drag.  The  last  term  in  Equation  (65)  acts  to  drive  the  angle  of  attack 
back  to  the  nominal  profile  (Harpold  and  Graves,  1978:125). 

The  control  algorithm  developed  above  is  implemented  subjected  to  a  roll 
constraint.  The  constraint  is  to  limit  the  bank  angle  rate  to  5  degrees  per  second.  Thus, 
each  time  the  control  algorithm  is  implemented,  a  commanded  bank  angle  is  determined. 
This  commanded  bank  angle  is  then  compared  to  the  current  bank  angle  to  see  if  it  is 
achievable  within  the  prescribed  constraint  of  5°/s  roll  rate. 

Roll  Reversal, 

As  the  vehicle  reenters  the  atmosphere,  the  bank  angle  will  be  modulated  to 
control  to  the  reference  drag  profile  and  to  null  range  errors.  In  order  to  reach  the  desired 
termination  point,  the  vehicle  will  have  to  reverse  the  direction  of  bank  at  certain  points 
along  the  trajectory.  The  prior  discussions  addressed  how  the  drag  acceleration  was 
affected  by  the  roll  reversal.  The  discussion  here  will  focus  on  when  roll  reversals  are 
implemented.  It  is  desirable  to  maintain  a  heading  angle  which  is  within  some  given 
deadband  of  the  desired  heading  angle.  Desired  heading  angle  is  calculated  as  the 
heading  angle  at  which  the  vehicle  would  travel  from  its  current  location  to  the 


49 


termination  point  along  a  great  eirele  route.  This  angle,  desired  azimuth,  AZd,  is 


ealeulated  from: 


AZj^  =  tan 


_ eos(Z^)sin(/^  -/J _ 

cos(4  )  sin(T^  )  -  sin(4  )  cos(T^  )  oos(/^  -  / J 


The  deadband  used  for  heading  angle  eontrol  is  as  shown  in  figure  26: 


(66) 
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Figure  26:  Azimuth  Error  Deadband 


The  actual  heading  angle  of  the  vehicle  is  calculated  from  Equation  (3)  and  compared 
with  the  deadband.  When  the  deadband  is  exceeded  a  roll  reversal  will  be  commanded. 
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IV,  Results 


Chapter  Overview 

This  chapter  will  present  the  results  obtained  when  applying  the  drag  profile 
guidanee  coneept  to  the  developed  SMV  model.  Two  mission  profiles  will  be  examined 
and  an  error  analysis  will  be  conducted. 

High-Inclination  Orhit 

The  first  case  considered  will  be  a  75°-inclination  orbit  with  no  crossrange 
requirement.  The  initial  orbit  altitude  is  500  km.  A  deorbit  burn  of  121  m/s  (397  fl/s) 
will  put  the  vehicle  on  a  Hohmann  transfer  ellipse  to  an  altitude  of  80  km,  with  reentry 
beginning  at  122.2  km.  The  initial  velocity  is  7.94  km/s  and  the  initial  flight  path  angle  at 
122.2  km  is  1.09°.  The  landing  site  is  at  Kennedy  Space  Center,  FL. 

Figure  27  shows  the  ground  trace  of  the  initial  75°  orbit,  along  with  the  location  at 
which  reentry  begins  at  122.2  km,  and  the  landing  site. 

Figure  28  is  a  plot  of  altitude  versus  distance  for  the  reentry.  During  the  initial 
stages  of  reentry,  approximately  2000  km  of  range,  the  plot  is  essentially  linear  as  the 
flight  path  angle,  y,  is  essentially  constant.  At  approximately  80  km  the  effects  of  the 
sensible  atmosphere  begin  to  slow  the  vehicle’s  descent. 
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Figure  27;  Ground  Trace  of  Orbit  (i=75°) 
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Figure  28;  Altitude  vs.  Range  (i=75°) 


The  increased  atmospheric  effects  at  about  80  km  altitude  can  also  be  seen  in  Figure  29, 
which  shows  the  variation  of  velocity  with  respect  to  altitude.  The  velocity  remains 
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nearly  constant  over  40  km  of  altitude.  Once  an  altitude  of  approximately  80  km  is 
reached,  the  velocity  begins  decreasing  in  a  near-linear  manner. 
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Figure  29;  Altitude  vs.  Velocity  (i=75°) 

Figure  30  shows  the  reentry  corridor,  the  reference  profile,  and  the  actual  trajectory  in 
the  D-V  plane. 


Figure  30;  SMV  Reentry  Profile  in  D-V  Plane  (i=75°) 
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During  the  initial  stages  of  reentry  the  heating  eonstraint  is  the  primary  eoneern. 
In  partieular,  the  region  of  eoneern  is  from  an  altitude  of  approximately  85  km  down  to 
an  altitude  of  75  km — this  eorresponds  with  the  region  at  the  beginning  of  the  eorridor 
where  the  aetual  drag  aeeeleration  level  is  signifieantly  higher  than  the  nominal  drag 
aeceleration  level.  At  approximately  6  km/s,  whieh  eorresponds  to  an  altitude  of  65  km, 
the  maximum  normal  load  beeomes  the  limiting  eonstraint. 

The  variation  of  the  bank  angle  with  respeet  to  veloeity  is  shown  in  Figure  3 1 . 


Figure  31;  Bank  Angle  vs.  Velocity  (i=75°) 

The  bank  angle,  (j),  is  modulated  to  control  the  drag  acceleration  level  and  to  null  range 
errors.  The  modulation  is  fairly  smooth  except  in  the  regions  where  the  guidance  shifts 
from  one  phase  to  the  next.  At  these  transition  points,  peaks  are  seen  in  the  bank  angle 
history  as  the  vehicle  attempts  to  reach  nominal  conditions  before  entering  the  next 
phase.  The  bank  angle  history  shows  that  the  vehicle  will  perform  4  roll  reversals  in 
order  to  stay  on  azimuth  for  the  termination  point. 


54 


Figure  32  shows  the  angle  of  attack  modulation  throughout  the  reentry,  compared 
with  the  nominal  angle  of  attack  profile.  The  commanded  angle  of  attack  agrees  closely 
with  the  reference  angle  of  attack  profile,  since  the  bank  angle  is  being  modulated  to 
control  drag  acceleration  levels  and  null  range  errors. 
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Figure  32;  Angle  of  Attack  vs.  Velocity  (i=75°) 


Figure  33  is  a  plot  of  the  actual  and  desired  azimuth  angle  during  reentry.  The 
actual  azimuth  angle  drifts  within  the  deadband  through  the  reentry.  Whenever  the 
deadband  is  exceeded,  a  roll  reversal  is  commanded  (as  shown  in  Figure  31). 
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Figure  33;  Azimuth  angle  vs.  Velocity  (i=75°) 


Figure  34  is  a  plot  of  the  flight  path  angle,  y,  versus  velocity.  The  flight  path 
angle  is  1 .09°  at  initial  reentry  and  decreases  to  a  minimum  of  0.06°  at  an  altitude  of  76 
km.  For  the  majority  of  the  reentry,  the  flight  path  angle  remains  between  0°  and  2°.  As 
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Figure  34:  Flight  Path  Angle  vs.  Velocity  (i=75°) 
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the  velocity  decreases  below  2  km/s  the  flight  path  angle  begins  decreasing  rapidly,  and 
the  vehicle’s  rate  of  descent  increases. 


The  primary  goal  of  the  guidance  algorithm  employed  in  this  study  is  to  guide  the 
vehicle  through  the  established  reentry  corridor  to  a  termination  point  at  an  acceptable 
energy  state  for  approach  and  landing.  Figure  30  shows  the  successful  guidance  through 
the  reentry  corridor  without  violation  of  constraints.  As  for  the  acceptable  energy  state  at 
the  termination  point,  Table  5  presents  a  comparison  of  the  actual  and  nominal  terminal 
values  for  velocity,  altitude,  and  position  for  this  scenario.  Figure  35  shows  the  ground 
trace  of  the  reentry  profile. 


Table  5:  Terminal  Conditions  for  SMV  Reentry  (i=75°) 


Conditions 

Nominal 

Actual 

Guidance  Error 

Velocity 

0.764  km/s 

0.764  km/s 

Okm/s 

Altitude 

29.4  km 

28.8  km 

-.6  km 

Latitude 

28.60°  N 

28.59°  N 

-.01° 

Longitude 

279.10°  E 

279.05°  E 

1 

o 

C/1 

0 

Thus,  the  vehicle  was  accurately  guided  to  within  3.8  km  of  the  termination  point  at  an 
altitude  of -0.6  km  below  the  nominal  desired  altitude,  without  violating  reentry 
constraints. 
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Figure  35;  Ground  Track  of  SMV  Reentry  Path  (i=75°) 

Low  Inclination  Orbit  with  Cross  Range  Requirement 

The  second  case  will  look  at  the  reentry  profile  for  a  vehicle  in  an  orbital  plane 
which  does  not  pass  directly  over  the  terminal  point.  The  initial  altitude  will  again  be 
500km,  thus  the  deorbit  bum,  initial  velocity,  and  initial  flight  path  angle  will  be  the  same 
as  calculated  in  the  first  case; 

AV=  121  m/s  =  397fps 
Vo  =  7.94  km/s  =  26000  ft/s 
Yo=  1.09° 

The  inclination  will  be  28.6°.  For  this  scenario,  a  hypothetical  orbit  at  28.6°  inclination 
which  passes  directly  overhead  of  the  terminal  point  will  be  referred  to  as  the  nominal 
orbit.  The  actual  initial  orbit  of  the  vehicle  will  be  rotated  22.5°  to  the  West  from  the 
nominal.  That  is,  the  right  ascension  of  the  ascending  node  of  the  actual  orbit  will  be 
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22.5°  West  of  the  right  ascension  of  the  ascending  node  of  the  nominal  orbit.  The  setup 
is  depicted  in  Figure  36.  This  scenario  is  significant  in  that,  for  a  rotating  earth,  it 
represents  the  single  orbit  (or  Abort  Once  Around)  scenario. 

Figure  37  shows  that  the  variation  of  altitude  with  respect  to  range  is  similar  to 
the  high- inclination  scenario. 


Figure  36;  Ground  Trace  of  Actual  and  Nominal  Orbits  (i=28.6°) 
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Figure  37;  Altitude  vs.  Range  (i=28.6°) 


In  Figure  38  it  can  be  seen  that  once  again  the  velocity  remains  essentially  constant  until 


an  altitude  of  approximately  80km. 
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Figure  38;  Altitude  vs.  Velocity  (i=28.6°) 


60 


Figure  39  shows  the  vehicle’s  reentry  profile  in  the  D-V  plane.  Once  again,  there  is  a 
period  of  increased  drag  early  in  the  profile,  but  the  heating  constraint  is  not  violated. 


Figure  39;  SMV  Reentry  Profile  in  D-V  Plane  (i=28.6°) 

Figure  40  shows  the  bank  angle  modulation  for  this  reentry  profile.  Since  the 
landing  site  does  not  lie  in  the  plane  of  the  orbit,  the  initial  azimuth  and  the  desired 
azimuth  are  not  the  same — as  shown  in  Figure  41. 
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Figure  40:  Rank  Angle  vs.  Velocity  (i=28.6°) 


Figure  41:  Azimuth  angle  vs.  Velocity  (i=28.6°) 


Figure  42  depicts  the  actual  and  nominal  angle  of  attack  profiles  for  this  scenario,  and 
Figure  43  shows  the  flight  path  angle  as  a  function  of  velocity. 
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Figure  42;  Angle  of  Attack  vs.  Velocity  (i=28.6°) 
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Figure  43;  Flight  Path  Angle  vs.  Velocity  (i=28.6°) 
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The  nominal  terminal  eonditions  for  the  energy  state  and  position  at  entry 
guidanee  termination  are  the  same  as  ealeulated  previously.  Figure  44  shows  the  ground 
traek  of  the  reentry  profile,  and  Table  6  summarizes  the  final  eonditions. 


Figure  44;  Ground  Traek  of  SMV  Reentry  Path  (i=28.6°) 


Table  6;  Terminal  Conditions  for  SMV  Reentry  (i=28.6°) 


Conditions 

Nominal 

Actual 

Guidance  Error 

Veloeity 

0.764  km/s 

0.765  km/s 

+.001  km/s 

Altitude 

29.4  km 

28.7  km 

-.7  km 

Latitude 

28.60°  N 

28.56°  N 

0 

P 

Longitude 

279.10°  E 

279.09°  E 

-.01° 

The  final  position  eorresponds  to  a  4.02  km  position  error  at  entry  guidanee  termination. 

Error  Analysis 

In  order  to  eomplete  the  analysis  of  the  applieation  of  the  developed  guidanee 
algorithm  to  the  SMV,  it  is  appropriate  to  perform  an  error  analysis.  Errors  in  the  initial 
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state  of  the  vehicle  will  cause  the  profile  to  be  shifted  off  the  nominal  profile.  A 
functional  guidance  system  should  be  capable  of  guiding  the  vehicle  to  the  termination 
point,  without  violating  entry  constraints,  in  the  presence  of  errors.  Possible  sources  of 
error  for  the  vehicle  initial  state  include;  altitude  at  guidance  initiation,  variations  in 
initial  orbital  parameters,  and  thruster  (deorbit  AV)  inaccuracies.  In  addition  to  errors  in 
the  initial  state,  variations  in  the  aerodynamic  properties  (Cl  and  Cd)  as  well  as 
atmospheric  variations  should  be  considered. 

For  the  purposes  of  error  analysis  in  this  study,  a  Monte  Carlo  analysis  will  be 
conducted.  “For  a  given  Monte  Carlo  trajectory  the  Gaussian  random  number  generator 
is  entered  with  the  appropriate  standard  deviation,  the  output  is  added  algebraically  to  the 
mean  value  and  the  trajectory  is  run”  (Regan,  1993:449). 

High  Inclination  Orbit  Error  Analysis. 

For  this  analysis,  the  profile  as  setup  in  scenario  1  above  will  be  used — high 
inclination  orbit  with  no  crossrange  requirement.  The  mean  values  for  use  in  the  Monte 
Carlo  simulation  will  be  the  nominal  value  for  each  initial  condition,  as  previously 
calculated.  The  standard  deviation  will  then  be  used  to  estimate  the  magnitude  of  the 
error  for  each  condition.  Table  7  summarizes  the  parameters  to  be  used  in  the  error 
analysis. 
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Table  7;  Initial  Conditions  and  Uncertainties  for  Monte  Carlo  Analysis 


Initial  Conditions 

Mean  Value 

Standard  Deviation 

Velocity 

7.94  km/s 

.01  km/s 

Flight  Path  Angle 

1.09° 

.01° 

Azimuth  Angle 

62.24° 

0.1° 

Altitude 

122.2  km 

1  km 

Latitude 

73.01°  N 

0.5° 

Longitude 

168.65°  W 

0.5° 

For  the  aerodynamic  coefficients,  the  mean  value  is  taken  to  be  the  value 
determined  from  the  curve  fits  in  Chapter  2  at  the  given  angle  of  attack,  and  the  standard 
deviation  is  taken  to  be  0.01.  For  atmospheric  variations,  the  standard  sea-level  density  is 
taken  as  the  mean,  and  the  standard  deviation  around  the  mean  is  0.07  kg/m^3  (Regan 
1993:449). 

Figures  45  and  46  show  the  results  of  a  100-run  Monte  Carlo  simulation  for  bank 
angle  and  angle  of  attack.  Figure  47  shows  the  reentry  corridor  in  the  D-V  plane  along 
with  the  profiles  followed  for  the  100  runs.  Some  variation  is  noted,  especially  near  the 
transition  points  from  one  guidance  phase  to  another,  but  the  reentry  constraints  are  not 
violated.  Figure  48  is  a  plot  of  the  termination  points  for  the  each  run.  The  mean 
distance  to  the  desired  termination  point  at  entry  guidance  termination  is  3.82  km,  with  a 
standard  deviation  of  1 .46  km.  The  maximum  distance  from  the  desired  termination 
point  is  10.79  km. 
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Figure  45:  Rank  Angle  vs.  Velocity  for  100-Run  Simulation 


Figure  46:  Angle  of  Attack  vs.  Velocity  for  100-Run  Simulation 
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Figure  47;  D-V  Plane  Reentry  Profile  for  100-Run  Simulation 
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Figure  48;  Entry  Guidanee  Termination  Loeation  for  100-Run  Simulation 
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The  distribution  of  the  altitude  at  termination  is  shown  in  Figure  49  with  the  mean 
altitude  at  entry  guidance  termination  being  28.98  km  and  a  standard  deviation  of  0.32 
km. 
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Figure  49;  Altitude  at  Entry  Guidance  Termination  for  100-Run  Simulation 
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V.  Conclusions  and  Recommendations 


Results  and  Recommendations 

The  objective  of  this  study  was  to  determine  a  means  of  guiding  a  conceptual 
Space  Maneuver  Vehicle  (SMV)  through  atmospheric  reentry  to  a  pre-determined  entry 
guidance  termination  point.  The  strategy  adopted  to  solve  this  problem  was  to  apply  the 
concepts  used  for  Space  Shuttle  Orbiter  reentry  guidance  to  the  SMV.  Having  done  this, 
the  following  conclusions  can  be  drawn  from  the  results  provided  in  Chapter  4. 

1 .  The  concept  of  selecting  a  desired  drag  acceleration  profile  and  controlling  the  vehicle 
to  that  profile  through  bank  angle  modulation  is  an  effective  method  for  guiding  the  SMV 
for  atmospheric  reentry.  The  results  for  the  two  scenarios,  as  well  as  the  100-run  Monte 
Carlo  analysis,  show  that  in  all  cases  the  vehicle  is  guided  to  the  termination  point,  within 
the  constraint  corridor.  In  all  cases,  the  energy  state  at  guidance  termination  compares 
favorably  with  the  desired  energy  state. 

2.  Although  the  entry  initiation  altitude  was  chosen  to  be  122.2  km,  atmospheric  effects 
are  essentially  negligible  until  the  vehicle  reaches  an  altitude  near  80  km.  Until  80  km 
altitude,  the  velocity  and  flight  path  angle  remain  essentially  constant. 
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3.  The  control  algorithm,  as  developed,  shows  some  weakness  in  the  region  near  the 
transition  points  from  one  guidance  phase  to  the  next.  The  discontinuities  at  the 
transition  points  in  the  nominal  drag  acceleration  profile  cause  the  controller  to  over¬ 
compensate  in  an  effort  to  match  range  and  drag  acceleration  levels  at  these  points.  As  a 
result,  “spikes”  in  the  bank  angle  and  drag  acceleration  profile  occur. 

4.  The  entry  guidance  algorithm  developed  is  designed  to  accommodate  prograde  orbits 
from  0°-89.99°  inclination.  Retrograde  orbits  were  not  studied,  however,  minor 
modification  to  the  algorithm  should  account  for  this  type  of  orbit. 

The  following  recommendations  for  further  study  are  given; 

1 .  The  analysis  conducted  herein  is  somewhat  limited  by  the  assumptions  made. 

Whereas  this  analysis  could  be  used  as  an  effective  baseline,  the  fidelity  of  the  simulation 
would  have  to  be  improved  for  an  operational  guidance  method.  In  particular,  a  more 
robust  model  of  the  Earth’s  atmosphere  accounting  for  cross-winds,  hemispheric 
variations,  and  seasonal  variations  should  be  input  into  the  guidance  algorithm.  Also, 

a  higher  order  model  of  the  Earth’s  gravitational  field,  to  include  variation  with  altitude, 
as  well  as  oblateness  effects  should  be  incorporated. 

2.  The  vehicle  parameters  used  herein  are  all  estimates  for  the  Space  Maneuver  Vehicle. 
As  the  vehicle  moves  from  the  drawing  board  to  the  prototype  phase,  these  parameters 
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will  become  more  obtainable.  The  parameters  of  mass,  wingspan,  and  area  may  not  be 
affected,  but  the  aerodynamics  are  likely  to  require  modification.  Construction  of  a 
higher- fidelity  model  than  the  one  used  herein  for  use  in  determining  the  aerodynamic 
coefficients  is  one  possible  approach.  Ideally,  as  a  prototype  is  developed,  wind  tunnel 
data  and  actual  flight  test  data  could  be  used  to  more  accurately  describe  these 
coefficients  in  the  entry  guidance  algorithm. 

3.  The  concept  of  controlling  the  vehicle  to  a  desired  drag  level  used  for  Shuttle  entry 
guidance  is  a  tried  and  tested  concept.  A  similar  approach  is  to  “define  the  reference 
profile  as  a  function  of  the  energy,  specifically,  a  piecewise-linear,  continuous  function  of 
the  energy”  (Lu,  1997:143).  Using  this  approach  preserves  the  analytical  range 
prediction  and  allows  for  non-computationally  intensive  on-board  determination  of  the 
optimal  reference  drag  profile  (Lu,  1997:143).  As  a  course  of  further  study,  this  concept 
could  be  applied  to  the  SMV  reentry  guidance  and  compared  with  the  results  obtained 
using  the  D-V  plane  entry  corridor. 

4.  Once  again,  the  goal  of  this  study  was  to  develop  an  entry  guidance  controller  for  the 
SMV — not  to  optimize  the  reentry  profile.  The  controller  developed  here  could  be  used 
to  determine  optimal  reentry  trajectories  in  a  further  study.  Studies  of  interest  would 
include  the  determination  of  optimal  profile  for  accumulated  heat  load  or  the 
determination  of  angle  of  attack  profile  to  provide  maximum  crossrange. 
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Appendix  A 


Measurements  are  scaled  to  the  vehicle  length:  L  =  29  ft  (8.84m) 


a.  0.1 74L 


b.  0.279L 


V 


c.  0.220L 


V 


d.  0.1 52L 


e.  0.175L 


Figure  50:  MATLAB  SMV  Model  Top  View  1 


a.  Nose  Length 

b.  Length  from  Nose  to  forward  lifting  surface 

c.  Forward  lifting  surface  inboard  length 

d.  Length  from  rear  of  forward  lifting  surface  to  front  of  rear  stabilizer 

e.  Length  from  front  of  rear  stabilizer  to  rear  of  vehicle 
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Measurements  are  scaled  to  the  vehicle  length:  L  =  29  ft  (8.84m) 


f  0.058L  g.  0.093L  h.  0.058L 


i.  0.238L 


◄ — ► 


◄ - ► 


◄ — ► 


<■ 


Figure  5 1 :  MATLAB  SMV  Model  Top  View  2 


f.  Rear  stabilizer  outboard  length 

g.  Rear  stabilizer  inboard  length 

h.  Forward  lifting  surface  outboard  length 

i.  Delta  wing  inboard  length 
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Measurements  are  scaled  to  the  vehicle  length;  L  =  29  ft  (8.84m) 


< - - ► 

k.  0.069L 

◄ - ► 

l.  0.163L 


-4  -6 


- >> 

m.  0.523L 


Figure  52;  MATLAB  SMV  Model  Front  View 


j.  Rotation  angle  for  Rear  Stabilizer 

k.  Nose  diameter 

l.  Body  width 

m.  Wingspan 
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Appendix  B 


Alternate  Ranging  Method 

An  alternate  method  to  caleulate  the  remaining  range  to  go  for  the  vehiele  is  to  extend  the 
integration  of  the  equations  of  motion  to  the  final  part  of  the  trajeetory.  For  this  portion 
of  flight  the  angle  of  attaek  is  set  to  20°,  which  closely  corresponds  to  the  maximum  L/D 
for  the  vehicle.  The  initial  flight  path  angle,  y,  is  set  to  15°,  and  the  initial  altitude  and 
velocity  are  29.4  km  and  0.764  km/s  (2500  fps),  respectively.  Range  is  maximized  by 
setting  the  bank  angle,  (j),  to  zero,  thereby  achieving  maximum  lift.  Integration  of  the 
equations  of  motion  subjected  to  the  above  constraints  yields  a  range  of  82.6  km.  The 
flight  profile  is  shown  in  Figure  53. 
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5 

oL - 1- - 1 - 1 ^ ^ ^ ^ - 1 

0  10  20  30  40  50  60  70  80  dO  100 

Ground  dstartce  (km) 


Figure  53;  Altitude  vs.  Distance  for  Terminal  Flight  ((j)=0°) 
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The  constraints  above  assumed  that  the  vehicle  would  fly  with  zero  bank  angle  to 
maximize  range.  This  will  likely  not  be  the  case  as  the  vehicle  maneuvers  into  position 
with  the  runway  for  landing.  To  estimate  the  range  capability  in  the  presence  of  non-zero 
bank  angle,  an  average  bank  angle  can  be  used.  This  approach  allows  an  estimate  of 
range  to  be  obtained  without  knowing  the  terminal  flight  path  in  advance.  Figure  54 
shows  the  effects  of  a  30°  average  bank  angle  on  the  range. 


Figure  54;  Altitude  vs.  Distance  for  Terminal  Flight  ((t)avg'=  30°) 
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Appendix  C 


Initial  Lat/Long  Determination 

Determining  the  latitude  and  longitude  of  the  initial  reentry  interface  point  is 
essentially  a  problem  in  spherical  trigonometry.  The  first  step  is  to  determine  the  reentry 
range,  R,  using  Equations  47-49  with  the  appropriate  initial  and  final  conditions  for 
velocity,  energy,  and  drag.  Next,  the  ground  trace  of  the  initial  orbit  is  determined  and 
plotted.  The  final  step  is  to  determine  the  point  in  the  orbit,  approaching  the  landing  site, 
that  lies  a  distance  R,  or  more  specifically  an  angle  0  =  R/Re,  from  the  landing  site.  This 
is  accomplished  by  choosing  a  longitude  at  least  180°  from  the  landing  site.  Then,  using 
the  accompanying  latitude  for  the  orbit,  the  range  angle  between  this  point  and  the 
landing  site  is  determined  using; 

R^ngie  -  *  [sin(T)  sin(T^  )  +  cos(T)  cos(T^  )  cos(/  -  /^  )]  (67) 

where; 

L  =  Latitude  from  ground  trace 
Ll  =  Latitude  of  landing  site 

I  =  Longitude  from  ground  trace 

II  =  Longitude  of  landing  site 

If  the  angle  is  greater  than  the  desired  range  angle,  then  the  longitude  is  incremented  and, 
using  the  accompanying  latitude,  the  range  angle  is  again  calculated.  This  process  is 
continued  until  the  range  angle  either  equals,  or  becomes  less  than,  the  desired  range 
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angle.  The  position  in  the  orbit  where  this  oeeurs  defines  the  latitude  and  longitude  of  the 
initial  reentry  interface. 
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